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Abstract. Spatially averaged inhomogeneous cosmologies in classical general 
relativity can be written in the form of effective Fricdmann equations with sources that 
include backreaction terms. In this paper we propose to describe these backreaction 
terms with the help of a homogeneous scalar field evolving in a potential; we call it the 
'morphon field'. This new field links classical inhomogeneous cosmologies to scalar field 
cosmologies, allowing to reinterpret, e.g., quintessence scenarios by routing the physical 
origin of the scalar field source to inhomogeneities in the Universe. We investigate a 
one-parameter family of scaling solutions to the backreaction problem. Subcases of 
these solutions (all without an assumed cosmological constant) include scale-dependent 
models with Friedmannian kinematics that can mimic the presence of a cosmological 
constant or a time-dependent cosmological term. We explicitly reconstruct the 
scalar field potential for the scaling solutions, and discuss those cases that provide 
a solution to the Dark Energy and coincidence problems. In this approach, Dark 
Energy emerges from morphon fields, a mechanism that can be understood through 
the proposed correspondence: the averaged cosmology is characterized by a weak decay 
(quintessence) or growth (phantom quintessence) of kinematical fluctuations, fed by 
'curvature energy' that is stored in the averaged 3-Ricci curvature. We find that the 
late-time trajectories of those models approach attractors that lie in the future of a 
state that is predicted by observational constraints. 



PACS numbers: 04.20.-q, 04.20.-Cv, 04.40.-b, 95.30.-k, 95.36. +x, 98.80.-Es, 98.80.-Jk 
1. Introduction 

The fact that the spatially averaged inhomogeneous Universe does not evolve as 
the standard model of cosmology, furnished by a homogeneous-isotropic solution of 
Einstein's laws of gravity, has recently become a major topic aiming at a possible 
solution to the Dark Energy problem [SU], (§], [Z31- While the averaging problem in 
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relativistic cosmology has a long history, initiated by George Ellis [33] (references may 
be found in [33] and JI]), the backreaction terms, i.e. those averaged contributions that 
lead to deviations from a Friedmannian cosmology, have been considered quantitatively 
unimportant. Recent work that conjectures a large backreaction effect (e.g. jUj) has 
been accompanied both by counter conjectures (e.g. [2]) and special solutions of the 
averaged Einstein equations [7] that support the claim that backreaction could be made 
responsible for the Dark Energy gap (e.g. [TU], [3JJ], [3J, [HID EH], [HD] and references and 
discussion in [HE]). 

In parallel to this discussion, other models have been extensively investigated, 
most notably quintessence models that invoke the presence of a scalar field source, 
this scalar field being a standard one or a phantom field (with a negative kinetic 
energy) [23 E21 EE1 IS] and references in the reviews 69] and (2HJ - Those models imply 
the possibility of a scenario featuring a time-dependent cosmological term. Particular 
properties of the scalar field potential are discussed on phenomenological grounds, which 
can lead to an acceleration of the regional Universe and so to an explanation of Dark 
Energy. A further problem in conjunction with modeling a repulsive component in 
universe models is known as the coincidence problem, i.e. a recent domination of this 
component is favoured, possibly occuring around the epoch of structure formation. Its 
solution also motivates the construction of models with a time-dependent cosmological 
term that needs no 'fine-tuning'. 

In this paper we also consider scalar field cosmologies, but we shift the perspective 
from the usual interpretation of a scalar field source in Einstein's equations to a mean 
field description of averaged inhomogeneities. We see a number of advantages entailed 
by such a description that motivated the present investigation. First, a (homogeneous) 
scalar field as a model of spatially averaged (i.e. effectively homogeneous) geometrical 
degrees of freedom that are physically present does not need to be justified as an 
additional source arising from fundamental field theories; it has a well-defined physical 
status and, as such, does not suffer from a phenomenological parameterization, since 
it is constrained by Einstein's field equations. Second, inhomogeneities encoded in 
backreaction terms give rise to the scalar field cosmology, fixing its potential, its equation 
of state, etc.; they not only influence the evolution of the scalar field, but they determine 
it. Third, the proposed correspondence allows a realistic reinterpretation of quintessence 
models and other phenomenological approaches involving scalar fields. These approaches 
must not be considered as independent alternatives, they here describe the same physics 
that also underlies the backreaction approach, and so can be confronted with physical 
constraints; in other words, if backreaction can indeed be made responsible for the late- 
time acceleration of the Universe, then the effort spent on quintessence models and other 
approaches [28J can be fruitfully exploited in terms of the proposed rephrasing. Fourth, 
since the underlying effective equations are general and do not invoke perturbative 
assumptions, the scalar field cosmology provides access to the whole solution space, 
notably to the non-perturbative regime. Finally, subcases of exact solutions include a 
cosmological constant, which therefore can or cannot be included in Einstein's equations; 
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the cosmological constant is a particular solution of the averaged inhomogeneous model 
without necessarily being present in Einstein's equations. 

Since the scalar field introduced in this way stems from averaged geometrical degrees 
of freedom, i.e. kinematical backreaction terms encoded in the extrinsic curvature of 
spatial hyper surf aces, as well as their scalar curvature, we call it the morphon field: 
it effectively models the inhomogeneities shaping the structure in the Universe and 
capturing the total effect of kinematical backreaction, i.e. the effect leading to a 
deviation from the kinematics of the standard model. Speaking in terms of geometry, we 
may roughly say that the morphon models the geometrical vacuum degrees of freedom 
in Einstein's theory, which are neglected when modeling the Universe with the help 
of averaged matter sources only. We call the scalar field a morphon only if it arises 
from the proposed mean field description; we can of course still consider other scalar 
fields as sources of Einstein's equations. We expect, and we shall demonstrate this 
for quintessence fields, that the morphon is capable of acting like any other scalar 
field model, e.g., assigned to an inflaton, a curvaton, a dilaton, etc., depending on our 
capability to exploit the proposed correspondence, i.e. the possibility to find solutions 
of the averaged Einstein equations and to reconstruct the scalar field potential for this 
solution. 

We here provide a first step that should help to appreciate the various possibilities, 
but this investigation does not claim to be exhaustive. In this line we are going to 
choose a minimal parameterization of the scalar field cosmology in the following sense: 
first, we do not include the cosmological constant and the constant curvature term 
in the 'Friedmannian part' of the averaged equations, since both arise as subcases of 
special solutions for the backreaction terms; second, we choose the simplest foliation 
(constant lapse and vanishing shift in the ADM (Arnowitt Deser Misner) formulation 
of Einstein's equations) and the simplest matter model 'irrotational dust' in order to 
exemplify the correspondence; third, we realize the correspondence for a particular 
class of scaling solutions with single-power laws. In this framework we show that the 
correspondence holds with a standard minimally coupled scalar field that can play the 
role of a quintessence field. If we include the full coordinate degrees of freedom as well 
as an inhomogeneous pressure source in a perfect fluid energy momentum tensor (or 
even imperfect fluid sources), more sophisticated scalar field cosmologies would arise. 
We comment on this possibility at the end of the paper. 

The proposed correspondence can also be inverted, e.g. we may start with a known 
model of quintessence and try to recover the corresponding solution of the averaged 
Einstein equations. We shall not exemplify this inversion in the present paper, we just 
note that free parameters in a given quintessence model can be determined through this 
correspondence. 

We proceed as follows. Section 2 recalls the basic equations and relations needed 
for the present work. Section 3 sets up the correspondence between kinematical 
backreaction and scalar field models. Section 4 investigates a family of scaling solutions 
to the backreaction problem. Section 5 exploits the proposed correspondence. Here 
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we explicitly reconstruct the potential for the scaling solutions, explore the solution 
space of inhomogeneous cosmologies with the help of the scaling solutions, and discuss 
some particular cases, which have been advanced before in the literature, and which 
are now interpreted with the help of a morphon field. We also give a concrete model 
and confront it with observational constraints. Section 6 contains a summary and an 
outlook on possible generalizations. 

2. The backreaction context 

The averaging problem in relativistic cosmology involves a variety of approaches and 
problems. Generally, research work in this field deals with averaging inhomogeneities 
in matter and/or geometry. Mostly, spatial averaging is envisaged, but also averaging 
on the light-cone is accessible by some averaging techniques (compare discussions and 
references in [35], [221 ECU EH ) 

In this work we focus on the comparatively simple approach of averaging the 
scalar parts of Einstein's equations on a given foliation of spacetime. The time- 
evolution of integral properties of the cosmological model on compact spatial domains 
can be extracted from Einstein's equation without any perturbative assumptions by 
Riemannian volume integration. The simplest example is the time-evolution of the 
volume. Thus, we do not aim at changing the physics of the inhomogeneous cosmological 
model. More ambitious and physically different in motivation are averaging strategies 
that effectively replace the inhomogeneous hypersurfaces and inhomogeneous tensor 
fields on these hypersurfaces by another, smoothed universe model [13], [211 ■ These 
latter techniques lead to further 'intrinsic backreaction' effects by, e.g. flowing averages 
on a bumpy geometry to averages on a constant-curvature 'template universe' as a 
fitting device. This can be nicely put into practice using the Ricci-Hamilton flow that 
renormalizes the averaged variables and leads to a 'dressing' of cosmological parameters 
by those additional 'backreaction' terms |lH].|14j. 

With kinematical averaging we aim at an effective description of the kinematics of 
the inhomogeneous Universe, and we still encounter a number of scalar contributions 
that add up to the averaged matter sources in an inhomogeneous model ('kinematical 
backreaction terms') as a result of the fact that spatial averaging and time-evolution 
are non-commuting operations (not necessarily to be attributed to the nonlinearity of 
Einstein's equations). Employing a perfect fluid source in the energy momentum tensor, 
those 'backreaction terms' comprise a contribution from averaged expansion and shear 
fluctuations (i.e., terms encoding extrinsic curvature of the hypersurfaces), the averaged 
3-Ricci curvature (i.e, a term that encodes intrinsic curvature of the hypersurfaces), 
an averaged pressure gradient (or an averaged acceleration divergence), and frame 
fluctuation terms, i.e. coordinate effects (like the averaged variance of the lapse function 
in the ADM setting). In this framework and for vanishing shift, but arbitrary lapse and 
arbitrary 3-metric, the general equations were given in [H|. However, we further restrict 
the analysis to a universe model filled with an irrotational fluid of dust matter [Zj in 
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order to provide the most transparent framework for the purpose of setting up the 
correspondence with a scalar field cosmology. 

2.1. Averaged ADM equations for constant lapse and vanishing shift 

We employ a foliation of spacetime into flow-orthogonal hypersurfaces with the 3-metric 
gij in the line-element ds 2 = —dt 2 + dX l dXK We define the following averager, 
restricting attention to scalar functions ^>(t,X l ): 

{^{t 1 X l )) v :=^- [ ^{t.X^Jd^X ; V v = [ Jd 3 X , (1) 
v v Jv Jv 

with J := ^det (gij); g^ is an arbitrary metric of the spatial hypersurfaces, and X 1 
are coordinates that are constant along flow lines, which are here spacetime geodesies. 
Defining a volume scale factor by the volume Vp of a simply-connected domain T> in a 
t-hyper surf ace, normalized by the volume of the initial domain T>\, 



the following exact equations can be derived [7j (an overdot denotes partial time- 
derivative). First, by averaging Raychaudhuri's equation we obtain: 

3— + 4irG(g) v -A=Q V ; (g) v = . (3) 

CLv VV-flv 

The first integral of the above equation is directly given by averaging the Hamiltonian 
constraint: 

2 



ouV SttG ! s A (TZ) V + Q V 



, (Q)v - T = - N /P J - , (4) 
a-p J 6 6 o 

where the total restmass Mp, the averaged spatial 3-Ricci scalar (1Z) V and the 
kinematical backreaction term Qx> are domain-dependent and, except the mass, time- 
dependent functions. The backreaction source term is given by 

Q v := 2 (II) V - ~ (1)1 = \{(e~ (6) v ?) v - 2 {a 2 ) v ; (5) 

here, I = Q t i and II = |[(0 l i) 2 — 0^-0^] denote the principal scalar invariants of the 
expansion tensor, defined as minus the extrinsic curvature tensor Kij := — 0y. In the 
second equality above it was split into kinematical invariants through the decomposition 
0y = \g%j9 + o~ij, with the rate of expansion 9 = Q\, the shear tensor a^, and the rate 
of shear a 2 := l/2aija^; note that vorticity is absent in the present foliation; we adopt 
the summation convention. 

The time-derivative of the averaged Hamiltonian constraint agrees with the 
averaged Raychaudhuri equation (j3J) by virtue of the following integrability condition: 

d t Q v + 6H V Q V + d t (K) v + 2H V (K) v = , (6) 
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where we have introduced a volume Hubble functional Hx> '■— af/ a v- The above 
equations can formally be recast into standard Friedmann equations for an effective 
perfect fluid energy momentum tensor with new effective sources jB]|: 

& = ^v-^Qv-^{n) v ; & = - 1 ± & Q v + ^L & (7l) v . (7) 

3— = A-4ttGG& + 3i&) ; 3H^ = A + 8nGgf s ; gf s + 3H V + pf s ) = . (8) 
a© 

Eqs. (JBJ) correspond to the equations @, (jlj) and (JUJ), respectively. 

Given an equation of state of the form pf s = j3(gf s ,aT>) that relates the effective 
sources (JZJ) with a possible explicit dependence on the volume scale factor, the effective 
Friedmann equations (jSJ) can be solved (one of the equations (JBJ) is redundant). 
Therefore, any question posed that is related to the evolution of scalar characteristics of 
inhomogeneous universe models may be 'reduced' to finding the cosmic state on a given 
spatial scale. Although formally similar to the situation in Friedmannian cosmology, 
here the equation of state depends on the details of the evolution of inhomogeneities. 
In general it describes non-equilibrium states. 

We finally wish to emphasize that these equations are limited to regular solutions: 
as in the non-averaged case, the matter model 'dust' generically leads to shell-crossing 
singularities. In the averaged equations this fact would be mirrored in a break of 
the boundary of the averaging domain or a merging of two boundaries (Legendrian 
singularities), thus inducing a jump of the Euler-characteristic of the boundary. This 
would especially happen for small collapsing domains and is related to the fragmentation 
and merging of structures. Here, by assumption, the domain T> must remain simply- 
connected. It is possible to cure this small-scale problem by generalizing the matter 
model; for example, if multi-streaming is accounted for, an extra term related to 
velocity dispersion would add up to the kinematical backreaction term. Such a term 
can regularize singularities as was discussed in detail within the Newtonian framework 
in [TSJ. 

2.2. Derived quantities 

For later convenience we introduce a set of dimensionless average characteristics in terms 
of which we shall express the solutions: 



v 8irG . v A v (K) v v Qv 

OXi-p dtlj) Dilp Dilp 



We shall, henceforth, call these characteristics 'parameters', but the reader should keep 
in mind that these are functionals on T>. Expressed through these parameters the 
averaged Hamiltonian constraint (JH) assumes the form of a 'cosmic quartet': 

n£ + nj + ng + nJJ = 1 . (10) 

\ Note that in this representation of the effective equations p e g denotes an 'effective pressure'; there 
is no pressure due to a matter source here. 
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In this set, the averaged scalar curvature parameter and the kinematical backreaction 
parameter are directly expressed through (71) v and Q v , respectively. In order to 
compare this pair of parameters with the 'Friedmannian curvature parameter' that is 
employed to interprete observational data, we can alternatively introduce the pair 



being related to the previous parameters by + figjv = fi^ + fig. 
These parameters arise from inserting the first integral of Eq. 

k v . 1 y* d 1 



(11) 



/ ^^^4(0 = ^(^ + 2^) , (12) 



into (gj): 



al + k Vi 8nG{g) v A_ 1 /•* . d 



°x> 3 3 3a|, j ; . 



/ dt' Q^4(t') . (13) 



This equation is formally equivalent to its Newtonian counterpart ^Bj, [IZ|- It 
shows that, by eliminating the averaged curvature, the whole history of the averaged 
kinematical fluctuations acts as a source of a generalized Friedmann equation. 

Like the volume scale factor a© and the volume Hubble functional H-p, we may 
introduce 'parameters' for higher derivatives of the volume scale factor, e.g. the volume 
deceleration functional 



^:=-^ir = ^ + 2fil-fi? , (14) 



or (volume) state finders [HHIE] (see also and references therein). 

In this paper we shall denote all the parameters evaluated at the initial time by the 
index V{, and at the present time by the index V . 

More details concerning these equations and their solutions may be found in [7| l%| IH| ITT] . 



3. The morphon field 

In the above-introduced framework we distinguish the averaged matter source and 
averaged sources due to geometrical inhomogeneities stemming from extrinsic and 
intrinsic curvature (backreaction terms). The averaged equations can be written as 
standard Friedmann equations that are sourced by both. Thus, we have the choice 
to consider the averaged model as a cosmology with matter source 'morphed' by a 
mean field that is generated by backreaction terms. We shall demonstrate below that 
this introduction of a 'morphon field' provides a natural description. We say 'natural', 
because the form of the effective sources in Eq. (JZJ) shows that, for vanishing averaged 
scalar curvature, backreaction obeys a stiff equation of state as suggested by the fluid 
analogy with a free scalar field jSH], [!]• Moreover, if we also model the averaged 
curvature by an effective scalar field potential, we find that the integrability condition 
© provides the evolution equation for the scalar field in this potential, and it is identical 
to the Klein-Gordon equation, as we explain now. 
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3.1. Setting-up the correspondence 

We now propose to model the effective sources arising from geometrical degrees of 
freedom due to backreaction (i.e., roughly the 'vacuum' part of the effective equations) 
by a scalar field $d evolving in an effective potential Ux> '■= £/($£>)§, both domain- 
dependent, as follows (recall that we have no matter pressure source here): 

Qes =■ {eh + q% ; pf s =■■ vl , (15) 

with [53], inn 

Ql = \&d + Uv ; P% = eU> 2 v - U v , (16) 

where e = +1 for a standard scalar field (with positive kinetic energy), and e = — 1 for a 
phantom scalar field (with negative kinetic energy). Thus, in view of Eq. (J7J), we obtain 
the following correspondence: 

1 Q v = e&l-Uv ; --^ (K) v = 3U V . (17) 



8ttG u ' 8ttG 

We see that the averaged scalar curvature directly represents the potential, whereas 
the kinematical backreaction term represents 'kinetic energy density' directly, if the 
averaged scalar curvature vanishes. This representation of <2x> is physically sensible, 
since it expresses a balance between 'kinetic energy' E£ n := 1/2$%,Vd and 'potential 
energy' E^ ot := -U V V V . For Q v = we obtain the 'virial condition' || 2eE^ in + E^ ot = 0, 
and so kinematical backreaction is identified as causing deviations from 'equilibrium' 
(defined through this balance^). Note that Friedmann-Lemaitre cosmologies correspond 
to the vanishing of <2x> (established 'virial equilibrium' or vanishing relative information 
entropy according to [13]) on all scales. The scale-dependent formulation allows to 
identify states in 'virial equilibrium' on particular spatial scales. Below we learn that 
this 'virial balance' can be stable or unstable. 

Inserting (JT7|) into the integrability condition (0) then implies that §x>, for &v 7^ 0, 
obeys the Klein-Gordon equation: 

^ v + 3H v ^ v + e^-U(^ T) ) = . (18) 
a<Px> 

With this correspondence the backreaction effect is formally equivalent to the dynamics 
of a homogeneous, minimally coupled scalar field. Given this correspondence we can 
try to reconstruct the potential in which the morphon field evolves. Note that there are 
two equations of state in this approach, one for the morphon, w$ := and the 

■P 



total 'cosmic equation of state' including the matter source term, wf s := p^/g^ 



§ We choose the letter U for the potential to avoid confusion with the volume functional. 

|| For negative potential energy (positive curvature) the sign e = +1, and for positive potentials the 

sign e = — 1 (phantom energy) is suggested from the scalar virial theorem 2E^ in = —E^ ot . 

% An alternative definition of 'out-of-equilibrium' states uses an information theoretical measure as 

proposed in |43j and discussed in the present context in 
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It is also worth noting that a usual scalar field source in a Friedmannian model, 
attributed e.g. to phantom quintessence that leads to acceleration, will violate the strong 
energy condition g + 3p > 0, i.e.: 

3- = -4irG(Q + 3p) = -4nG(g H + Qz + 3p*) >0 , (19) 
a 

and actually also the weak energy condition g + p > 0, while for a morphon field both 
are not violated for the true content of the Universe, that is ordinary dust matter. It is 
interesting that we can write a 'strong energy condition' for the effective sources, i.e.: 

3 ^ = _4ttG(^ + 3p? ff ) = -4ttG ( {g) v + g^ + 3p v )= -AnG(g)v + Q v < . (20) 
a v 

While we do not need 'exotic matter', the above condition will be 'violated' in order to 
have volume acceleration, Qv > AnG{g)v jUj, [H], cf. Subsection 15.4. II 



3.2. Newtonian limit 



In the Newtonian limit ^H], the above correspondence persists. The sources of the 



morphon in the effective Friedmann equations 
Eqs. (HU) and JEJ): 



are then identified as follows, cf. 



87TG 



V 
P<i>N 



1 



24vrG 



2 {t' 
d 



3k 



•V 



l a v 



(21) 

However, Newtonian cosmologies suppress the morphon degrees of freedom on some 
fixed large scale where the kinematical backreaction term has to vanish identically [16J. 
In particular, this remark applies to cosmological N-body simulations: by construction, 
these simulations enforce 'virial equilibrium' of the morphon energies on the scale of the 
simulation box. 

The Newtonian framework also offers a concise explanation of our choice 'morphon 
field': the kinematical backreaction term Q-p can be entirely expressed through 
Minkowski Functionals of the boundary of the averaging domain |9,. These 
functionals form a complete basis in the space of (Minkowski-) additive measures for 
the morphometry of spatial sets. 



3.3. Motivation: the morphon modeling a cosmological constant 

Let us give a simple example. As shown in jTT], Subsect. 3.2 (also advanced as a 
motivation case in jUj), the effective source Qt> may act as a cosmological constant. 
The general condition for the corresponding exact solution of the effective Friedmann 
equations © (with (g) v = and A = 0) is (R) v = Qk v .Ja v — 3Qx>', Qv = Qvv e.g., 
for kv { = we have: 

$© = ; g% = U v ; pI = -U v ; U v = U Vi . (22) 
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The kinematical backreaction term assumes a constant value Qz> = 8TcGUx> i} and the 
morphon potential mimics a (scale-dependent) cosmological constant. Note that the 
averaged curvature is non-zero, (TZ) V = —IAkGUt)-^ so that simultaneously the morphon 
unavoidably installs a non-zero averaged scalar curvature. 

We shall now exemplify this correspondence for a new class of scaling solutions that 
contains this and other known subcases. 



4. A family of scaling solutions of the backreaction problem 

In the following, we shall study a class of solutions that prescribes backreaction and 
averaged curvature functionals in the form of scaling laws of the volume scale factor a-p. 
If not explicitly stated otherwise, we restrict attention to the case A = throughout, 
and treat the cosmological constant as a particular morphon. 



4-1. Exact scaling solutions 

In this subsection we shall present a systematic classification of scaling behaviors for 
the cosmological models introduced previously. The averaged dust matter density (q) v 
evolves, for a restmass preserving domain V, as (q)d = (f?)r>i &t>~ Let us suppose that 
the backreaction term and the averaged curvature also obey scaling laws, that is: 

Qv = Qv, <%> ; (R) v = K Vi a p v , (23) 

where Qz>. and TZd, denote the initial values of Qx> and (R) V) respectively. 
Rewriting the integrability condition (JHJ), 

vd t {a v Q v ) + a v 2 d t (4 (K) v ) = , (24) 



a 



a first scaling solution of that equation is obviously provided by ([7j Appendix B): 

Qv = Qv, a v & ; (K) v = K Vl a^ 2 . (25) 
Moreover, this is clearly the only solution with n ^ p. In the case n = p, we define 

a new 'backreaction parameter' r v (that can be chosen differently for a chosen domain 

of averaging + ) such that Q©. oc 7£z>.; the solution reads: 

Qv = r (K) v = rH v ,a v , (26) 

where (cf. Fig. |TJ) 

n = -2£±M , (27) 
(1 + r) 

with r 7^ —1. The case r = — 1 is not represented as a solution in this class; this 
line of states degenerates to a point corresponding to a model with Einstein-de Sitter 
kinematics, i.e. it has vanishing backreaction and vanishing averaged curvature. Note 
here, that the vanishing of Qt>, if required on all domains, is necessary, but also sufficient 
for the reduction of the averaged equations to Friedmann-Lemaitre cosmologies. This 
limiting case also appears in the relations among the cosmological parameters, discussed 
after the next subsection. 



+ 



For notational ease we henceforth drop the index T> and simply write 



r. 
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n 




-2 -1,5 -1 -0,5 0,5 1 

r 



Figure 1. Scaling law n for the backreaction as a function of the parameter r. 
4-2. Discussion of the scaling solutions 

The solution (J25|) corresponds to the case where the backreaction and the averaged scalar 
curvature evolve independently, leading to an averaged curvature similar to a constant 
'Friedmannian curvature' and an additional term that scales as (f?)|>. As the universe 
model expands, this solution is rapidly equivalent to a pure constant curvature term, 
we may say that this solution represents (and maintains) a near-Friedmannian state 
(Qv decays much more rapidly oc V^ 2 compared with the averaged density oc V^ 1 ). 
Therefore, we shall not model this solution by a scalar field when describing the late- 
time dust-dominated Universe. It is interesting to note that this solution exhibits the 
same late-time behavior as the long-wavelength part of the solution found through the 
gradient expansion approximation scheme in |47j . 

On the contrary, the solutions entail a strong coupling between kinematical 
backreaction and averaged scalar curvature. The coupling itself must be considered a 
generic property; it has been identified as being responsible for a much slower decay of 
kinematical fluctuations in an expanding universe model on the cost of averaged scalar 
curvature. This is a genuinly relativistic property. It was argued in ^0] that it is this 
possibility which is needed for an explanation of Dark Energy through backreaction 
today, provided the initial conditions are appropriate. 

The scaling solutions (pjj) can be employed to represent generic features of 
backreaction- or curvature-dominated cosmologies (while density fluctuations must 
not be large). Such cosmologies may significantly deviate from a standard Friedmann 
cosmology with regard to the temporal evolution of their parameters. We note that, even 
if deviations in the volume scale factor ax> at a final time may not be large, deviations in 
in its history, i.e. its time-derivatives, in particular its second time-derivative related to 
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the cosmological parameters, may be large (this insight is a result of a detailed analysis 
in Newtonian cosmology [TTj). 

Let us make a general remark concerning the solution subspace of the scaling 
solutions given above. We appreciate that the polynomial nature of (j2SJ), together with 
the form of the integrability condition (JHJ), implies that any linear combination of these 
solutions provides a new solution. In particular, one can always add a constant curvature 
term to a particular solution. We also infer that only the case r = (n = —2) implies 
a (scale-dependent) 'Friedmannian' evolution of the (physical) curvature parameter 
Q^. Only the singular case r = — 1 would imply a vanishing 'Friedmannian' curvature 
parameter (see below). If we would require r = — 1 on all scales, then the model reduces 
to the Euclidean case (everywhere vanishing 3-Ricci curvature). 

We resume this discussion in the next section with the help of concrete examples. 
There we also provide an illustration and further discussion of the subspace defined by 
the scaling solutions. 

Finally let us note that, while the scaling solutions and scenarios investigated 
here and below satisfy the averaged equations of motions, these models are still 
phenomenological and not derived non-perturbatively from the fundamental theory, i.e. 
there is no guarantee that corresponding realistic solutions of the original inhomogeneous 
Einstein equations could be found that satisfy the assumed scaling laws. 

4-3. Some relations among cosmological parameters 

We here write some useful relations among the dimensionless cosmological parameters, 
as they were introduced earlier, Eqs. (jHI) and ((TTJ) . For the scaling solutions (J27|) we 
have: 



and for the 'Friedmannian curvature parameter' we find (by integration of (j!2)l for the 
solutions (1271) ): 



(The latter equation follows by noting that = 0,^ {H^J H"^) a£> and 4r/(l + r) = 
— (n + 2)). Eq. ()29|) now explicitly shows that the generally held view that Q]? models 
the averaged curvature is a misperception: as soon as kinematical backreaction is 
relevant (itself or its time-history), the averaged curvature may evolve very differently 
compared with a constant-curvature model. Evaluating (|29|) at initial time, we find 
that initial data differ only by the parameter r in this class of solutions, which eases 
their observational determination: 



Q _ 






(29) 




(30) 
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Since in the sum + fig = &q N + fij^, the 'Newtonian' parameter VLq Ni cf. Eq. (fTTj) . 
would be interpreted as a cosmological constant parameter, fijj', in a 'Friedmannian 
fitting model', we can already from these relations infer that with = —kvj (H^a^), 
fi^ 1 = (1 + r)fi^, and fif = fi^ 1 (H^, /H^a^), the value of a fitted A— parameter would 
directly depend on the initial datum for fi^ according to the relation 

+ = + = VLf^al = fi£ a™ +2 , (31) 

so that today 

ng> + fig = of + fil? = a n+2 (32) 

With this relation we can determine the dependence of the scaling solution parameter 
r on the cosmological parameters today, expressed through the 'Friedmannian fitting 
parameters' fij^ and fij^ . This allows us to put constraints on r, which will be done in 
Subsections 15.4.21 and 15.4.31 

We see that in these estimates we also need the correct evolution of a-p in order to 
calculate au . This is possible for the particular class of scaling solutions we consider 
(we have numerically evaluated the evolution for the volume scale factor below); in 
general it requires a detailed model for the evolution of inhomogeneities. 

Finally, we emphasize that the above relations are a result of our 'single-scaling' 
model ansatz. More flexible models arise by superimposing scaling solutions (and so 
modeling the dynamics more realistically). The relation ()31|) is a consequence of this: 
in general, initial data for the 3-Ricci curvature and the integration constant /cx> ; can 
be independently chosen. The reader should therefore make up their mind about a 
particular choice of superposition of scaling solutions that they would like to implement. 

In this line our model is maximally conservative concerning the amount of 'early' 
Dark Energy [3T], i.e. the value of fig^y, interpreted as a (constant) A— parameter at 
the present time evolves as a time-dependent 'cosmological term'. Its initial value can 
be calculated: 

^qn = + fig " = (1 + r)ti% - ^ = . (33) 

Thus, we implicitly require the initial contribution of this term to vanish, while actually 
a value in the range of a few percent would be allowed by observational constraints [36J, 

HE], d- 

5. Morphon— quintessence 

The obvious candidates for scalar field models, which come into the fore in the situation 
of a dust-dominated Universe, are quintessence models. We shall now exploit the 
proposed correspondence to explicitly reconstruct the scalar field dynamics. Since 
quintessence models aim at mimicking a repulsive component in the cosmological 
evolution, it is natural that a 'working model' of quintessence would, via the 
proposed correspondence of a 'morphon-quintessence', also lead to a 'working model' of 
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backreaction. For the purpose of concretizing the correspondence, let us now consider 
solutions of the type ()26|) . We will show that they can be put in correspondence with 
a one-parameter family of homogeneous scalar field solutions that act as standard or 
phantom quintessence fields. 

5.1. Reconstructing the potential of the morphon field 

The scaling solutions (ffilj) provide: 

Kb, (_ , 1\ . . „ Kb, , _ „(l + 3r) 



This correspondence defines a scalar field evolving in a positive potential, if 7£x>. < 
(and in a negative potential if IZ^ > 0), and a real scalar field, if e7^x>,(r + 1/3) < 0. In 
other words, if TZ-d-, < we have a priori a phantom field for r < —1/3 and a standard 
scalar field for r > —1/3; if TZx> l > 0, we have a standard scalar field for r < —1/3 and 
a phantom field for r > —1/3. 

The system (|34|) can be inverted to reconstruct the potential of the morphon field. 
Indeed, the kinetic term of the scalar field can be expressed in terms of d^x>/dax>, and 
this equation can be explicitly integrated, for (g)vi 7^ 0, leading to: 



Moo ) = Vgj±gg+g arsinh f ./-(»+;>** 




l-3r)v^G 167rG(e)2?, 

V / 0^h^) > (35) 



arsinh 



(n + 3)V7rG 

where we have defined the fraction 7^ m of the curvature and density parameters: 

V ._ ^1 _ X>i (n+3) . / I q\ — (1 ~ 3r ) . B, ._ /o fi X 
•— — lnm a V i \n + 0)— ^ +r ^ ) Iflm ■— ^T>i ' 

In this relation, necessarily, r / -1/3 and r 7^ —1. We immediately find that <&t>{o>v) 
is an increasing function of a©. Then, inverting that relation and inserting the result 
into the expression for the potential, we obtain the explicit form of the self-interaction 
term of the scalar field: 



1 + r)K Vi f„ . , ,* ^-4<i±£> / (1 - 3r)V7rG 



= 3i+I> ((1 + r)o£,)^ -hrfftl (<^§V2^,) , (37) 

where (q)v\ is the initial averaged restmass density of dust matter, and where the 
restrictions introduced above still hold, with the new constraint r ^ 1/3. To sum up, in 
order to obtain a consistent description in terms of a real-valued scalar field, we must 
have TZv i > and e = +1 for r < —1, and TZx> i < for r > —1, with e = —1, if 
— 1 < r < —1/3 and e = +l, ifr>— 1/3. One can immediately notice that the energy 
scale of the scalar field potential is determined by the averaged matter density: the 
scales that determine the scalar field dynamics are fixed by the matter distribution. As 
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a result of our correspondence to quintessence models and in view of many results that 
were obtained in this field, the above potential can also be found in [SH] (their Eq. (121) 
with a typo corrected in [SB!*) an d> e -S- F^! - 

We now study the cases that were excluded in the above derivations. First, let us 
consider the case of vanishing matter source. The correspondence given by equations 
(J33|) and ()37|) holds in the presence of a non-vanishing matter field, but one can also 
reconstruct the scalar field cosmology for the vacuum. Setting (f?)©. = in the equations 
and (fT3*|) . and applying the same procedure as the one described above, one finds: 



of matter, this is exactly the solution (|3*Tj) in the case $x> — > +oo. 

As noted above, the morphon field can be interpreted as representing the effect 
of the averaged geometrical degrees of freedom; then the comparison of its potential 
in the presence of matter and in the vacuum tells us, how the matter field influences 
the backreaction terms: it affects the energy scale through a simple factor depending 
on the initial averaged matter density, and so modifies the dynamics of the morphon 
when the domain volume is small (because $£> is an increasing function of a©; the limit 
$73 — > +oo corresponds to a© — ■> +oo.) On the contrary, when the domain volume 
becomes big, the dynamics of backreaction is similar to its dynamics in vacuum, which 
is natural because the averaged matter density is then diluted. 

Second, we discuss the three cases of the backreaction parameter r that were not 
considered until now: r = ±1/3 and r = —1. In the case r = 1/3, the solution reads 
Qx> = 1/3 (R) v oc a^ 3 ; this corresponds to a scale-dependent Einstein-de Sitter scenario 
with a renormalized initial dust density (§)■£>.— 1Zt>J3GiiG (c/. Appendix A). This model 
has zero effective pressure p e s. The case r = —1/3 leads to Qx> = — 1/3 (H) v = const., 
which is equivalent to a scale-dependent Friedmanian scenario with a cosmological 
constant: A = Q©. = —TZvjS (that was our motivating example). The case r = — 1 
corresponds to a strict compensation between the kinematical backreaction and the 
averaged scalar curvature. It leads to a scale-dependent Friedmann model with only a 
dust matter source, in other words, to a scale-dependent Einstein-de Sitter model. 

The above three cases appear as limiting cases of the scalar field model. 

The scaling solutions correspond to specific scalar field models with a constant 
partition of energy between the kinetic and the potential energies of the scalar field. 
Indeed, if we define the kinetic energy by E^ in := |<&|,Ve> and the potential energy by 
-^pot := —U($t>)Vv as before, we find the following 'balance condition' for the scalar 
field representation of backreaction and averaged scalar curvature in the case of the 

* Thanks to Varun Sahni, who has pointed this out to us. 




(38) 



Up to the renormalization factor 
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scaling solutions: 

EL + {1 ^El t = . (39) 



We previously discussed the case r = ('zero backreaction') for which this condition 
agrees with the standard scalar virial theorem. This balance between kinetic and 
potential energies is well known in the context of scaling solutions of quintessence (see 
|5()| E2j and references therein) . 

Finally, the effective equation of state for this morphon field is constant and given by: 
„ l(l-3r) 1, 

»S = -37TTO = -3< B + 3 > ' (40) 

which is less than — 1/3, iff r e] — 1; 0[. The overall 'cosmic equation of state' including 
the matter source term is given by: 

W CS V ~ O l^r - W$ _ (n+3) V , l 4i J 

Q <* 3 (1 + r)a^ + l/ 7 ^ 1 + < VTfcm 



Or, equivalently: 



with 7& := = (1 + r)7£„ . (42) 



W ? ff = <(1-^) . (43) 

Thus, the 'cosmic state' asymptotically evolves, for an expanding universe model with 
(n + 3) = (1 — 3r)/(l + r) > 0, into wf s — > — {n + 3)/3 = A necessary condition 
for the scalar field part to dominate the expansion of the universe model at late times 
is n = —2(1 + 3r)/(l + r) > —3, which implies r < 1/3. In that case, since the 
equation of state w$ is constant, a universe model dominated by backreaction and 
averaged curvature approaches the following evolution of the volume scale factor (note 
civ-, = a v (ti) = 1): 



a v (t) = (l + C(t-t,)) (1+r)/(1+3r) = (l + C(t-t t ))- 2/n ; 

with C : = i^±^^-(l + r )7^/6 = -yv^ • (44) 

For the class of scaling solutions we can obtain the time-evolution of the volume scale 
factor, and we appreciate that the scaling solutions explore possibilities similar to a 
Friedmannian evolution with a cosmological constant. Figure 121 illustrates this for some 
chosen values, where the density parameter is held fixed for comparison with our later 
analysis. For this figure we numerically integrate Eq. (£Q) for A = and with the scaling 
solutions inserted: 

(a v \ 2 8nG(g) Vi (1 + r) n 

W " "3""? 6~ {n)via " 

8ttG (q) Vi 2 / , 23, n n 

~ ~3 ofT" ' v ~ Vi V"o|" fc D 
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Figure 2. Volume scale factor ar> as a function of t for various scaling solutions; note 
that all the models shown have vanishing cosmological constant. The models had been 
integrated with a Vo = 1000, H = 70 km/s/Mpc, fl%> = 0.27 and Qq° + fl^ = 0.73. 
A classification scheme comprising the various cases will be provided below. Within 
this scheme the first row displays examples of Case D, Case C and the Einstein-de 
Sitter solution; the second row an example of Case B for a model of quintessence, the 
solution mimicking a cosmological constant, and an example of Case A for a model of 
phantom quintessence. 



The integration is performed using a non-stiff predictor-corrector Adams method 
provided by Scilab, with a Vo = 1000, H Vo = 70 km/s/Mpc, Q%° = 0.27, and 
fig° + f^° = 0.73, which yields for the constants in the previous equation: 8irG(g)'p i /3 = 
al o H Vo nl° and (1 + t)TZ v J& = -a^H Vo (Q% + flg>). 

5.2. Classification of the scaling solutions 

The following classification summarizes the results: 

• Case A: for r g] — 1; — 1/3[ we have a phantom scalar field with a positive potential 
of the form Ux> oc sinh /3 (a$x>) with (3 > 0; these models may lead to an accelerated 
expansion. A particular example of this type is analyzed below. 

• for r = — 1/3 we have a model that corresponds to a scale-dependent cosmological 
constant given by A = Qd. = —1/3TZt> ; - 

• for r > —1/3 we have a standard scalar field. Here, we can distinguish various 
cases: 

— Case B: for r g] — 1/3; 0[, the potential is positive and of the form sinh ai ($£>), 
with aii g] — 4;0[, leading to a quintessence field that may produce an 
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accelerated expansion. At the beginning, when $£> ~ 0, the potential 
is an inverse power-law one, corresponding to the so-called Ratra-Peebles 
potential [HZJ U(® v ) oc $ v a2 , with a 2 = 4(1 + 3r)/(l - 3r); when $ v 
becomes big, the potential is equivalent to an exponential potential, also well- 
known in the quintessence context |S3 EI]: U(<&v) oc exp(— a^v), with 
a 3 = 2^/(1 + 3r) /(l + r). It should be noted that the problem emphasized 
in jSUJ E2] for this kind of potential (that is the necessarily small scalar field 
density because of primordial nucleosynthesis constraints) doesn't hold here, 
since this potential arises during the matter-dominated result of 

backreaction, so that in such a scenario the scalar field is significantly sourced 
only during late stages of the matter-dominated era. The more r approaches 
— 1/3, the more this field mimics a cosmological constant. 

— Case C. for r G [0;l/3[ the potential still behaves as sinh Q4 ($£>), but with 
a 4 < —4. This potential is too stiff and the model cannot produce an 
accelerated expansion (w$ > —1/3). 

— for r = 1/3, the model is equivalent to a standard Einstein-de Sitter model 
with a scale-dependent and renormalized initial dust density and zero effective 
pressure p e s (cf. Appendix A). 

— Case D: for r > 1/3, the potential is of the form sinh a5 ($£>) with a 5 > 4 and 
the model cannot produce an accelerated expansion. 

• Case E: for r < —1, we have a standard scalar field rolling in a negative potential 
that is not bounded from below. Whereas such scalar fields are pathological 
when considered like fundamental scalar fields, this solution may be physical in 
the backreaction context. Indeed, the four preceeding models all correspond to 
7^x>i < 0, and this one corresponds to TZ^ > 0. We expect that more realistic 
solutions, modelled e.g. by a suitable superposition of scaling solutions, could 
provide potentials with minima, hosting 'bound states'. 

5.3. The solution space explored by the morphon 

The various cases listed above appear as 'cosmic states' that separate 'cosmic phases', 
illustrated in the following phase diagram, Fig. El To understand this diagram, we 
remark that the solution space for the effective cosmologies that are sourced by dust 
matter and 'morphed' by the scaling solutions form one-dimensional subsets in the two- 
dimensional space that is defined by Hamilton's constraint (taken at fixed spatial scale 
V and at A = 0), + VLq + = 1- Also, (scale-dependent) 'Friedmannian' models, 
characterized by the backreaction parameter r = and A = 0, form a one-dimensional 
subset defined by Vt^ + = 1. With the scaling solutions we are also restricted to 
measure zero sets, but we have a one-parameter family of them which allows us to 
explore the solution space. 

To illustrate the solution plane, we plot, instead of fig (related to Q^) the volume 
deceleration parameter q v = 1/2Q^ + 2fig, Eq. (JHJ), as a function of the only free 
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parameter fl^ in the 'cosmic phase diagram' Figure El For the scaling solutions we 
simply have: 

q l + r 2(1 +r) m 2 2 m 1 ' 

A priori, the scaling solutions can describe the whole plane, including cosmologies 
with > 1. Nevertheless, > 1 implies, from Hamilton's constraint |T0|) with 
A = 0, that (1 + r^lZ-cn > 0> which is exactly the opposite condition to the one that 
holds in the correspondence with a real-valued scalar field. In other words, a real- 
valued morphon field is only defined for < 1, and we shall concentrate on this class 
of models in the following. 

Since VLq + would be interpreted, in a Friedmannian 'fitting model', as 
+ i w ith negligible k-parameter in the concordance model [4*81 . we can infer the 



6 A ' k 

corresponding value for a fitted A— parameter in the same diagram, since for negligible 
A;— parameter, VLq + = 1 — ~ Q^, where with the upper index F we refer to a 
'fitted' A— parameter. 

Another feature in the following diagram are arrows that illustrate the time- 
evolution of the respective parameters. Again, the models corresponding to > 1 are 
not analyzed here (among them are 'Big-Crunch-models', cf. Eq ()47)1 : their dynamics 
depends strongly on initial conditions). There are attractors and repellorsjj in this 
diagram. Most notably, for the cases of interest to us that feature a late-time 
acceleration, 'Friedmannian' states are repellors, i.e. 'near-Friedmannian' states evolve 
away down to the attractor solution fl^ = 0, where backreaction- (or curvature-) 
domination is completed. 

Let us specifically discuss the various cases. We denote values in the solution plane 
by (q v ; Vt^ n ) and concentrate in the following only on expanding universe models. We 
write: 



2 1 + c 2 a£+ 3 ' m " l + c 2 a™ +3 



qv(av) = - — — ^ ra+3 ; V m = - - , (47) 



with ci = 4rc2/(l + r), c 2 = (1 — and ax> being an increasing function of 

time. These formulae are helpful to determine the evolution of a particular solution 
in the phase plane (g^;^). We learn from Figure 01 that no model corresponding 
to Case C can produce acceleration. According to Eqs. (|H|), scaling solutions in the 
sectors corresponding to negative potentials (case E) and pit-type potentials (case D) 
are attracted towards the Einstein-de Sitter model (1/2; 1). In all the other cases (A, 
B and C), the Einstein-de Sitter model appears as a repellor, and the attractors are 
located on the line (q v = 2r/(l + r); fl^ — 0) for r G] — 1; 0]. Each point of the straight 
line r = 1/3, that corresponds to a renormalized Einstein-de Sitter scenario, is a fixed 
point. 

[j Compare the analyses in and |72) . Due to the proposed correspondence to a scalar field cosmology 
we can directly use the results on scaling properties of the scalar field investigated by J25J and 0]. In this 
context, the determination of equations of state from similarity symmetries also provides an interesting 
tool |74j . A full-scale investigation of a dynamical systems analysis is not provided here. 
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Figure 3. This 'cosmic phase diagram' is valid for all times and on all scales, i.e. 
it can be read as a diagram for the corresponding parameters 'today' on the scale of 
the observable Universe. All the scaling solutions are represented by straight lines 
passing through the Einstein dc Sitter model in the center of the diagram (1/2; 1). 
The vertical line corresponding to (q v ; 1) is not associated with a solution of the 
backreaction problem; it degenerates to the Einstein-de Sitter model (1/2; 1). This 
line forms a 'mirror': inside the cone (Case E) there are solutions with > 1 that 
cannot be related to any real-valued scalar field, but are still of physical interest in the 
backreaction context (models with positive averaged scalar curvature). Models with 
'Friedmannian kinematics', but with renormalized parameters form the line r = 1/3. 
The line r = are models with no backreaction on which the parameter fi]? varies 
(scale-dependent 'Friedmannian models'). Introducing A would just shift the whole 
diagram down. Below the line r = in the 'quintessence phase' we find effective models 
with subdominant shear fluctuations (Qx> positive, fig negative). The line r = —1/3 
mimics a 'Friedmannian model' with scale-dependent cosmological constant. The line 
below r = —1/3 in the 'phantom quintessence phase' represents the solution inferred 
from SNLS data (c/. Subsection l5.4.3|l . and the point at (q-^fi^) = (— 1.03;0) locates 
the late-time attractor associated with this solution. 
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Hence, the Einstein-de Sitter model is a saddle point for the scaling dynamics and 
small inhomogeneities with Qx> > should make the system evolve away from it. The 
sign of Qv is important: for all the models corresponding to r > or r < —1, that is the 
cases C,D and E, which cannot produce accelerated expansion, we have Qx> < 0. In other 
words, the kinematical backreaction is dominated by shear fluctuations, cf. Eq. (|HJ). 
This does not necessarily mean that the universe model is regionally (on the scale T>) 
anisotropic, because in these cases kinematical fluctuations decay strongly. On the other 
hand, cases A and B that could be responsible for an accelerated expansion correspond 
to Qx> > and have subdominant shear fluctuations. Therefore, these models can be 
regionally almost isotropic, although kinematical fluctuations have strong influence. 

Finally, recall that we have TZx> { < for r > — 1 and TZx> i > for r < — 1. 

5-4- Construction of a realistic model: estimation of parameters and initial conditions 

The scalar field behavior (i.e. the late-time behavior of the cosmological model in 
the cases r < 1/3) essentially depends on the value of the backreaction parameter r 
that describes the ratio between kinematical backreaction and averaged curvature (or, 
in the scalar field language, the ratio between the kinetic energy of the field and its 
potential energy). We are going to estimate this ratio from observational constraints on 
the equation of state. Before we do so, let us summarize the conditions relevant for a 
late-time behavior featuring 'volume acceleration', i.e. a-p > 0. 

5.4-1- Acceleration conditions The condition for an accelerating patch T> (which we 
are going to take as large as our observable Universe) follows from the averaged 
Raychaudhuri equation © |Ej, [TTj : 

Q v > AixG(q) v , (48) 

implying with (JHJ) and ((Hj): 

-«S > x ; f < o , (49) 

and, for the class of scaling solutions (|2T|) (fi^ > 0, fl^ = 0, and the definition (j3T)jO: 

ng=rng ,i.e., -r 7 £ m > \ ^ <f = ^Q + 2r 7 £ m ) < , (50) 

This condition is met, if (now inserting the solution (}36|) and the relation to the 
Friedmannian curvature parameter (J42|0 : 

V r Vi (n+3) _ 1 /r-n 

-ri nm = -j^-^lk^ a v > 4 ■ ( 51 ) 

A realistic model would meet this condition at some time in the evolution leading 
thereafter to the observed acceleration value, i.e. ideally at the epoch around the time 
of structure formation solving the coincidence problem. 

From what has been said above, such realistic cases require r < —1/3, i.e., n > 0; 
in the limiting case r = — 1/3 (n = 0) (which is the solution mimicking a cosmological 
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constant) the condition (joTj) reads 7^,a|> > 1/2, and with = ti^' 1 a v 2 and 
K = *v we find: 

r = ~ , n = : 7 L4 > \ , i-e. , t£ > ^ ■ ( 52 ) 

The marginal case n = —3 in the exponent of the volume scale factor (|51|) plays a 
particular role that we shall explain more in detail in Appendix A. In the context of 
the acceleration conditions we can gain a better understanding of the cases n > — 3 by 
the following remark. In this marginal case the averaged density and the kinematical 
backreaction have identical decay rates with respect to the volume scale factor in an 
expanding universe model. This means that, in order to get a positive acceleration at 
the present time, already the initial data must satisfy the conditions (pfK|) and ({HI)) . The 
necessary initial value for kinematical backreaction is then large and suggests that we 
are looking at a region D that is close to satisfy the condition needed for the stationarity 
of the cosmos. Note, however, that the strict single-scaling solution n = —3 (r = 1/3) 
does not admit acceleration, cf. Eq. (jlS} , but a superimposed regional fluctuation would 
admit acceleration (or deceleration) (cf. Appendix A). 

In the cases n > —3, the averaged density decays faster than kinematical 
backreaction; hence, to attain sufficient acceleration today, the model needs the less 
magnitude of kinematical backreaction the weaker its decay given in terms of n. Let 
us take a case close to n = — 2 (the case n = — 2 degenerates to r = 0, cf. Eq. f|2Tjl ). 
then kinematical backreaction can be three orders of magnitude weaker initially, if the 
scale factor advanced to a value of av = 1000 today. This remark makes clear that 
the solution sector n > —2 (r < 0) contains solutions that can potentially explain 
the Dark Energy problem even when starting with small expansion fluctuations at the 
CMB (Cosmic Microwave Background) epoch. Moreover, it contains solutions which 
also solve the coincidence problem, although a more natural solution would not be an 
exact scaling solution, but one that would inject more backreaction at the formation 
epoch of structure. However, for models with fl^ < 1, cf. Eq. (HHJ), we have to go to 
values of n ~ (r ~ —1/3) in order to find sufficient acceleration. 

It should be emphasized that the interesting sector n > —2 is not what we could 
find in a weakly perturbed FRW (Friedmann Robertson Walker) model. These states 
rely on a strong coupling of kinematical fluctuations to the averaged scalar curvature 
of the universe model. Kinematical backreaction can only decay at such weak rates (or 
even grow for r < —1/3), if the time-evolution of the averaged scalar curvature largely 
deviates from the time-evolution of a constant curvature model; intuitively speaking, 
averaged fluctuations are strongly supplied by the 'curvature energy reservoir'. 

5.4-2. Observational constraints on parameters Recall that the envisaged class of 
single-scaling solutions implies that our parameter choices are unambiguous: we only 
have to specify an initial condition, say or f^ 1 , the backreaction parameter r, and 
the value of the volume scale factor today a-p . 
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The following estimates are done on the assumption that the effective model is 
based on solutions for a matter density source and a morphon field that is realized by 
the particular class of scaling solutions discussed in this paper. We focus on the effect 
of the morphon field for vanishing cosmological constant, and would like to demonstrate 
that an observer who is using a Friedmannian 'template universe model' would interprete 
this effect by a cosmological constant today. Thus, we are forcing the effective evolution 
of the volume scale factor to match with that of a Friedmannian model with A at initial 
and final time. However, all of our parameters are evolving according to the 'best-fit' 
scaling solution in the averaged inhomogeneous model. In particular, this implies that 
the time-derivatives of the volume scale factor evolve very differently compared with 
a standard Friedmannian model. Thanks to the existing constraints on the standard 
Friedmannian models, as for example Cold Dark Matter models with a Dark Energy 
component that has a constant equation of state, this procedure reduces the number of 
parameters that we have to estimate. Indeed, we then only have to determine r and one 
value for the initial data, e.g. Q-^. 

We emphasize that the interpretation of observational data and the resulting 
constraints strongly depend on model assumptions, i.e. it is commonly assumed that 
a standard ACDM model is the correct one. This implies, in particular, that we are 
not constraining the parameters of the averaged inhomogeneous model by observations 
reinterpreted within the inhomogeneous cosmology. For example, since the spatial 
curvature is not constant, the formulae for angular diameter and luminosity distances 
cannot be taken as the FRW ones. Furthermore, we point out that this reinterpretation 
is challenging, for all the other observational predictions that are based on a perturbative 
approach, like large-scale structure characteristics, must be reconsidered. It is not 
obvious, and this work together with others (e.g. jHE] and references therein) provides 
plausible counter arguments, that the late Universe could be described by a perturbed 
FRW model, even if smoothed over large scales. 

In this sense, our analysis is a demonstration of what the observed values of the 
standard cosmological parameters would imply for the averaged quantities. 

With these assumptions the observer with a 'Friedmannian template' then faces 
the following relation today: 

fi^ -|- f^Q = f2y^ -|- , 

where the latter corresponds to the 'biased' interpretation of the true dynamics. The 
parameter r is fully specified by the energy content of the Universe today, since by 
Hamilton's constraint 0%,° + = 1 ~ fi m°- 

Directly following from the relation (p?Tj) we have: 



with (n + 2) = —± ^ ; 7A° k --=7± ■ (53) 

n + 6 In a Vo ilK Vt^° 

Note that, again by Hamilton's constraint, the condition 1 + 7^ > implies 

(l-^°)/^° > 0, which allows for two positive curvature today negative) 
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Figure 4. The backreaction parameter r is shown as a function of fiP° and fij[, 
as the effect of backreaction would be interpreted in a Friedmannian 'template', with 
av ~ 1000. 

with > 1, or a negative curvature today with < 1. For our purpose of fitting 
the inhomogeneous model to a Friedmannian 'template' with cosmological constant, we 
choose the latter option. Furthermore, for > 0, there is still the possibility that 
< — This last condition is clearly not satisfied in the late-time Universe, so in 
the following, we will restrict the analysis to > 0. 

Giving initial data from WMAP [7TJ, fi£ ~ 0.72 and ~ 0.28, and taking 
a Vo ~ 1000, Eq. determines the value of n to be close to but slightly larger than 
0, i.e. r slightly below r = —1/3 pointing to a phantom quintessence. However, taking 
into account a variation of the initial data, we detect a large sensitivity of the precise 
value for r to the initial data in our scaling solutions. This is shown in Figure EJ We 
therefore employ an orthogonal observational constraint on the Dark Energy component 
below. 

Figure 0] shows r as a function of and f2^ for the bounds given by WMAP 
third year data jllj: e]0; 0.006] and Vt F h = 0.72 ± 0.04, and for a Vo = 1000 that 
roughly corresponds to a range of solutions integrated from approximately the epoch of 
the CMB (Cosmic Microwave Background). 

We infer that when tends to 0, r tends to —1, and to when tends to 
+oo. Moreover, we notice that, while increasing from zero, the morphon field is 
first a phantom field, then a standard quintessence field. Finally, r is almost insensitive 
to Q^, and it depends strongly on only in the case when is very small. (Note 
that the latter is a 'pathology' of our single-scaling ansatz; this behavior could be cured 
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by superimposing a constant-curvature solution.) 

5-4-3. Constraining r by supernovae observations The preceding considerations giving 
r in terms of and (or in terms of and by replacing with 1 — Q^fQ^ 
in Eq. (JHHJl) can be used to determine r. We now estimate r through the Dark Energy 
equation of state that best matches supernovae observations. Then, we use relation 
to compute the resulting value of fi^ , in order to check the consistency with our fitting 
model. 

There is a simple way to estimate the backreaction parameter r, by requiring that 
the equation of state for Dark Energy cf. Eq. (jlQ)) , matches the one inferred from 
supernovae observations, denoted by wsNia, provided this one is constant. Then: 

_ 1 1 + SwsNla 
3 1 - WsNla 

Taking into account the last SNLS best fit for a flat Friedmann model sourced by Dark 
Energy with a constant equation of state WsNia ~ —1.02, = 0.73 and = 0.27 
jS], we find r ~ —0.34, again suggesting a phantom quintessence. This value does 
not depend sensitively on variations in WsNia, and is therefore a more robust estimate 
compared with our previous one. 

Once this ratio is fixed, we can find the values of the initial data. We assume that, 
today, + fig = 0^ and find, e.g. for the initial ratio of the curvature parameter to 
the density parameter: 

7 ^ = WnG(g) Vi = (l + r )^a& o +3) ' 

which, at the CMB epoch, approximately setting a Vo to the value ax> ~ 1000, is 
7 ^~2.710- 9 . 

Here, we can check the consistency of the SNLS fitting with the curvature 
assumption. Indeed, in the SNLS fitting, we assumed that = 0, and, using the 
inferred value for r, r ~ —0.34 in Equation (J53)) . we find that Vt^° ~ 510~ 7 f2j[, in 
accordance with the assumption = 0. This shows that r as determined only through 
the Dark Energy equation of state is compatible without further assumptions. 

If we define an 'effective redshift' through the volume scale factor as in Friedmann 
cosmology, we can derive the effective redshift at which the expansion accelerates. It 
corresponds to a 'cosmic equation of state' for matter plus backreaction wf s ~ —1/3. 
Inserting ()41j) into this relation, we find an acceleration scale factor and an effective 
acceleration redshift: 

a& t =(-l/{4 n 8j)^ - 569 =► = — -1-0.76. 

\ / ax> 

The scalar field behavior and the potential corresponding to this model, as well as the 
time-evolution of various cosmological parameters, are presented in Figs. El and 
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Figure 5. Upper left: scalar field $p as a function of the volume scale factor a-p, 
in units of 2\J-kG. Upper right: Potential U(<&v), in GeV m~ 3 . Central left: volume 
scale factor as a function of time in Gyr. Central right: deceleration parameter as a 
function of time in Gyr. Lower left: the 'cosmic equation of state' w e s as a function of 
time in Gyr. Lower right: 'Friedmannian curvature parameter' as a function of time 
in Gyr, featuring a maximum at the onset of acceleration. All solutions are calculated 
for the SNLS best fit model. 
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Figure 6. The dimensionless cosmological parameters are shown as a function of time. 
While decays from 1 to the present value ~ 0.27, the sum ~t~ ^72, 

grows from 

to the value ~ 0.73, interpreted as the cosmological constant parameter today, fi£, in 
a Friedmannian 'template model'. Kinematical backreaction Qx> is positive (dominant 
expansion fluctuations) and grows slightly, implying a negative backreaction parameter 
that becomes more negative. This growth is supplied by the negative averaged 
curvature that evolves strongly, i.e. the positive curvature parameter grows from 
almost to a value larger than 1 at the present epoch. This behavior of the 'physical' 
curvature parameter should be compared with the evolution of the 'Friedmannian 
curvature parameter' displayed in the last figure. All solutions are calculated for the 
SNLS best fit model. In this model the parameters (f2^; J7g; 51^) evolve towards 
their attractor (0; —0, 515; +1, 515), where backreaction- and curvature-domination is 
completed. Note that the starting values differ from (Q^ 1 = l;f2g' = 0; Q^' = 0); if 
they would not, then the model evolves without backreaction as the standard model, in 
which the latter initial data set is maintained in the course of evolution. The standard 
model is a repellor for the scenario shown. 
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6. Concluding Remarks and Outlook 

Our proposal of a mean field description of backreaction effects through a minimally 
coupled 'morphon field', does not only provide a rephrasing of the kinematics of 
backreaction in terms of a scalar field cosmology, but it also justifies existence of the 
latter due to the fact that we identify averaged inhomogeneities in Einstein gravity as 
the underlying fundamental physics. Thus far, such a well-defined link was missing; 
alternatively, it is thought and there exist some plausible hints that the low-energy 
scale in some contemporary particle physics models would provide this link (e.g. [30J, 
[SSI, ESI), i-e. that the scalar field emerges from the dynamical degrees of freedom 
stemming from extra dimensions ('moduli fields'). 

The set of spatially averaged equations together with the mean field description of 
kinematical backreaction by a morphon field shows, in particular, that the framework 
of Friedmann-type equations is very robust. We understand now that an effective and 
scale-dependent Friedmannian framework is applicable with the surprising new input 
that, besides averaged matter sources, also a scalar field component emerges and has 
not to be invented. 

6.1. Summary 

We have exploited the proposed correspondence by investigating a family of scaling 
solutions of the backreaction terms. The study of scaling solutions is well-advanced 
in research work on scalar field models, and our study allows to translate those 
results into the backreaction context. Here, our discussion was not exhaustive. We 
concentrated on exploring the solution space of inhomogeneous cosmologies with the 
help of scaling solutions with a 'minimal' parameterization. Therefore, as a next 
step, this parameterization could be expanded and congruences with other work on 
scalar field models could be worked out. Note that an obvious such expansion would 
analyze superimposed scaling solutions that could be modeled either again by a single 
effective scalar field, or else by multiple scalar fields [27|, [22] • One example of a 
superposition of scaling solutions is provided in Appendix A. Furthermore, observational 
constraints set on quintessence or phantom quintessence models have direct relevance to 
constraints that have to be set on morphon fields [H|. However, observations have to be 
reinterpreted within the inhomogeneous cosmology underlying the morphon dynamics. 

In making the correspondence concrete we have not touched the question whether 
realistic dynamical models of the inhomogeneous Universe would comply with the class 
of averaged models that we singled out as 'realistic case studies'. For the scaling solutions 
the scalar field correspondence revealed, that the constancy of the fraction of kinetic 
to potential energy of the morphon field relates to averaged models that are driven by 
strong coupling between averaged 3-Ricci curvature and kinematical fluctuations, the 
constancy implying their direct proportionality. While in the scalar field context this 
assumption is often employed, here we get an interesting picture of what this represents. 
The key for a physical interpretation of the proposed correspondence is the equivalence 
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of the integrability condition (0) and the Klein-Gordon equation (JTSJ): 

a v 6 d t (a 6 v Q v ) +a v 2 d t (al(TZ} v ) = , & $ v + 3H v $v + e^-U(^ v ) = . 

Although the resulting picture may appear complex, it can be structured by looking 
at a single player: the averaged 3-Ricci curvature. If it is positive, it represents an 
energy reservoir, stored in a negative potential in the correspondence. If this reservoir 
exceeds the 'virial energy', then the excess energy is converted (in the scaling solutions 
directly) into an excess of kinetic energy, i.e. kinematical backreaction. The averaged 
curvature then decays faster than the 'Friedmannian' curvature. On the other hand, 
if it is negative, then the same logic applies: the averaged curvature consequently 
becomes stronger negative compared with the 'Friedmannian' curvature. Solutions of 
the Dark Energy problem, where a large positive value of backreaction is needed, will 
have to make strong use of this conversion, while 'near-Friedmannian' models don't. A 
Newtonian or post-Newtonian approach suppresses these degrees of freedom by freezing 
the averaged curvature to the 'Friedmannian' value (compare an attempt to work with 
the Friedmannian curvature parameter that, of course, needs an extra scalar field as a 
source for Dark Energy [HH]). From these remarks we understand the importance of a 
strongly evolving averaged 3-Ricci curvature, if the attempt to explain Dark Energy 
through backreaction should be successful. In this line there is clear support for the 
importance of a strongly evolving 3-Ricci curvature from the Lemaitre-Tolman-Bondi 
solution, see [03], and in particular the recent work [60J and references to other papers 
on this solution therein. Recently, Rasanen [00] provided an illustrative example and a 
comprehensive discussion of the physics of backreaction-driven accelerated expansion. 

Weakly perturbed Friedmann models have a special status if placed into the full 
solution space of the averaged inhomogeneous models: only if the averaged fluctuations 
decouple from the averaged curvature, i.e. if they evolve independently, then the 
averaged curvature evolves as a constant 'Friedmannian' curvature, and fluctuations 
decay in proportion to the square of the inverse volume. Hints that point to the 
likely existence of a curvature-fluctuation-coupling come again from averaging the LTB 
solution [HH] (see also [2S]) and other global solutions [TT] that all show an extra term in 
the averaged scalar curvature. In these examples the extra term evolves in proportion 
to the inverse volume, hence deviates from a 'Friedmannian' curvature and implies 
maintainance of large kinematical fluctuations that decay only in proportion to the 
inverse volume. While this particular behavior of the backreaction terms still requires 
large backreaction to start with ^0], e.g. a feature of globally stationary solutions 
[TT] . a stronger evolution of averaged curvature, i.e. injecting more curvature energy 
into kinematical fluctuations, would allow to start with 'near-Friedmannian' initial 
conditions and still explain current observations. An example for the latter possibility 
has been analyzed in this paper. This possibility can be interpreted such that there is 
some hope to find enough backreaction by starting with almost FRW initial data, as 
suggested by |37| . although here the conversion of curvature energy into backreaction 
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must be very efficient; a more moderate evolution into backreaction-dominated phases 
would result, if we already start with 'out-of-equilibrium' initial data. 

In this line we have also pointed out, and this is worth stressing again, that the 
'Friedmannian' curvature parameter is in general unrelated to the physical averaged 
scalar curvature, it is an integration constant and actually an integral of motion [11J: 
the dynamics of the averaged physical curvature is not represented by and, in most 
examples, deviates substantially from it. The restriction to a 'Friedmannian' evolution 
of the averaged curvature singles out the special case where the curvature-fluctuation- 
coupling is strictly absent. 

Summarizing these thoughts we can say that there are two conceivable scenarios 
that remove the need for Dark Energy: first, a 'soft scenario' that was discussed in detail 
in ^l] . Here we already have an initial global state with strong expansion fluctuations, 
so that regionally a moderate evolution of backreaction and averaged curvature suffices 
to explain the observed acceleration. It was, however, pointed out that such models 
imply paradigmatic changes, and observational data have to be reinterpreted in order 
to put firm constraints on Q-p at the CMB epoch. Second, there is a 'hard scenario' 
(implied by the suggestion of Kolb et al. jUj) that literally creates enough backreaction 
out of 'nothing', cf. Eq. (|33|). In our example of a particular scaling solution a phantom 
quintessence scenario arises that complies with the strong energy condition and also with 
constraints in accord with those already put on the standard model. Consequently, this 
scenario needs strong evolution of backreaction and averaged curvature as seen best for 
the dimensionless cosmological parameters in Fig. |H1 The present work has shown that 
the 'hard' version can be made consistent with the framework of the averaged Einstein 
equations. 

6.2. Outlook 

There are obviously a number of possible routes for generalizing the scalar field 
correspondence. Let us briefly discuss some of them. 

Spatial averages are scale-dependent, since we integrate over a given compact 
volume of the space sections. In this paper we left the domain-dependence untouched, 
all of our considerations were focussing on a fixed scale. However, the scalar field 
correspondence holds on every scale, which is not only reminiscent of, but also physically 
a manifestation of renormalizable quantities. In this respect the morphon analogues of 
quintessence are different from standard quintessence models. In this context we could 
'Ricci flow' the averages to averages on a constant curvature geometry [13], or we could 
study other renormalization group methods to control the scale-dependence, e.g. [20J, 

021, EH- 

A possibly fruitful investigation would consider string-motivated gravitational 
theories and, employing the proposed correspondence, would aim at determining the 
scalar field theory from the higher-dimensional geometry of extrinsic and intrinsic 
curvature. For this end one would have to derive the averaged equations for the 
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extended spacetime. Brane world cosmologies could be analyzed in the spirit of this 
correspondence. Not only in this context it will be interesting to understand, in which 
cases a non-minimally coupled morphon field would arise. 

Furthermore, there are other interesting strategies of a more technical nature that, 
however, widen the physical context of applications. For example, the introduction 
of a complex scalar field, which is necessary to access solutions with > 1 with a 
positive averaged 3-Ricci curvature; the discussion of globally stationary cosmologies 
in shows that a stationary cosmos obeys these conditions at early times, and later 
the averaged scalar curvature decreases fast until it becomes negative, while maintaining 
large kinematical backreaction. Another example is the averaging on a different foliation 
of spacetime, i.e. introducing coordinate degrees of freedom. The latter is actually 
necessary if matter sources other than 'irrotational dust' are considered; it results in an 
extra coupling of the scalar field to the matter sources stemming from a backreaction 
term due to the pressure gradient that is absent in the present paper |S]. This is the 
subject of a follow-up work that will also allow access to realistic models of backreaction- 
driven inflation Here, the present investigation in principle allows to model 

inflationary scenarios too by translating typical characteristics of inflaton fields into 
corresponding morphon fields, however, within the restricted cases of 'dust matter' or 
'vacuum'. 

Since $p models the inhomogeneous 'vacuum part' of the sources, there might be 
an interesting connection to the energy budget of gravitational waves that could be 
fruitfully exploited [25] . and, since <3>x> represents a 'mean field' of fluctuations, also a 
connection to statistical thermodynamics is implicit. In order to establish this latter 
link, however, we have to note that Qx> models the averaged spatial variance of extrinsic 
curvature and so far not 'fluctuations' in the thermodynamical sense. Another intimate, 
but on the level of a 'dust matter model' formal, relation to effective thermodynamic 
models is suggested in the case of imperfect fluid models. They imply non-equilibrium 
effects that in turn can be associated with a 'friction coefficient' proportional to the 
'cosmic equation of state' p c s/ Q c s- Thermodynamic arguments that were used in [£2|I7T)] 
within the imperfect fluid picture lead to further possibilities of interpreting the different 
scaling regimes, e.g. by employing the second law of thermodynamics the authors of 
|H2J ED! would conclude p e s < 0, and thus Qt> > 1/3(7Z)t>, i.e. r > 1/3 for positive 
averaged scalar curvature and r < 1/3 for negative one. Since the first case does not 
give rise to acceleration one would be led to conclude that the scalar curvature has to be 
negative on average, a conclusion shared by other work (in Appendix A we also reach 
this conclusion and provide references). 

A final comment related to scalar field perturbations, often investigated in a post- 
Newtonian setting, is in order. In those perturbative approaches, the scalar degree of 
freedom that emerges is a combination of matter and metric inhomogeneities. Within 
the well-developed standard perturbation theory ]§T\, it is natural to first restrict the 
analysis to a tight range around a FRW background in order to calculate backreaction 
effects PE], [HUE]. The amplitude of scalar perturbations, but also their derivatives 
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and, hence, all curvature terms have to be small jH]. Besides considerations of long- 
wavelength perturbations (e.g., , the analysis of sub-horizon perturbations with 
the aim to explain Dark Energy in a perturbative setting is the focus of many recent 
research papers and, here, we do not enter into the details of these works. We wish to 
point out two aspects that emerged from the present investigation, and which may be 
relevant for future strategies in a perturbation analysis. 

First, we address the averaging issue: a post-Newtonian approximation may (and 
in most cases of cosmological relevance will) be adequate locally or piecewise on a small 
range of scales. However, as soon as we are looking at integral properties, i.e. integrating 
out a wider range of scales, its applicability should be verified. In standard perturbation 
theory spatial averages are taken with respect to a 'background observer'. Since a major 
player in the mechanism that can produce large kinematical fluctuations is the averaged 
3-Ricci curvature, we have to be careful in relating Euclidean 'background averages' to 
the Riemannian volume averages that govern the dynamics of the averaged cosmology. 

Second, in the present approach we do not specify a metric of the space sections; 
the formulation and the correspondence holds for arbitrary 3-metrics. Also, when we 
spoke about a 'scale-dependent Friedmannian model', we referred to the kinematics of 
the volume scale factor, we did not refer to the FRW metric. We can, however, specify 
a spatial metric to establish a dynamical model. Here, we think that the metric setting 
must allow for large deviations of the 3-Ricci curvature from the constant Friedmannian 
curvature. As a 'rule of thumb' (another was recently given in [SI] for super-Hubble 
perturbations), we may say: if the averaged scalar curvature evolves at or near the 
constant curvature model, then there is no hope to model cases that lead to enough 
late-time acceleration. 

A future strategy related to perturbation theory could be motivated by Newtonian 
cosmological models for structure formation: in the Newtonian framework, an Eulerian 
perturbation theory does not provide access to the highly non-linear regime. Instead, 
the Lagrangian point of view offers a way to move with the largely perturbed fluid. 
A relativistic Lagrangian perturbative approach is currently worked out aiming at 
generalizing the Newtonian work ^7] that has employed the exact averaged equations to 
construct a non-perturbative model for inhomogeneities out of perturbatively calculated 
fluctuations. A Lagrangian perturbation scheme itself does not include non-perturbative 
features, that may be needed in this context; non-perturbative effects have been recently 
discussed in the Newtonian framework 12 . 

All these efforts aim at constructing a generic evolution model. The morphon field 
is an effective description without any perturbative assumption. But, one would like to 
establish the underlying inhomogeneous dynamics in order to understand, to what the 
realistic case studies that 'explain away' the Dark Energy problem actually correspond. 
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Appendix A: a relation to Priedmannian fitting models and the globally 
stationary solution 

A particular scaling behavior of the solutions 1)27)1. in which kinematical backreaction 
and averaged scalar curvature are proportional, can be exploited to define a mapping 
of the solutions ()2T)l for the particular case r = 1/3, i.e. n = —3, to a Friedmannian 
cosmology. 

Before we explain this mapping let us recall that the case r = 1/3 also arises for 
the stationarity condition, required in JT] for the global scale (extending the domain T> 
to the whole (compact) manifold E, and setting A = 0): 

3 as = _4vrG(,) s + Q Si = Q 

This condition implies either a globally static model or a globally stationary model 
featuring the solution: 

Qfe - -3- , (/<-) 2 - -2 + —3- , (A.3) 

where = C/a^ with: — 6C 2 = IZ^ — 3Qs ; . 

The relation to a Friedmannian 'fitting model' arises by noting that, in order to 
obtain the above solution, we do not need to assume the stationarity condition ()A.2)1 . 
Indeed, inserting the scaling solutions (|27)l into the averaged equations (J3j) and (JIJ) on 
any given domain T>, and superimposing a constant curvature solution to the averaged 
curvature, we obtain (restricting again attention to A = 0): 

a v , ^G(g) Vi -Q Vi n ( a v V I^G{q) v . - (J + l)Q Vi C 



3— + x 3' Ui = ; ^ (\ 3 — - — = • (A.4) 

We find that, for the special case r = 1/3, we can simply redefine the initial constants, 
fc£:=-C 2 ; 47cG(q}v. ■= 4:1tG(q} t >. - Q v . ; + 0£ := + ^ + fig = 1 , (A.5) 
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where f2^ F and f2jP denote the resulting dimensionless functionals which we are going 
to use as Friedmannian 'fitting parameters', so that Eqs. (jA.4|) assume the form of a 
constant-(negative) curvature Friedmannian model. 

Let us exemplify this correspondence. Setting ~ in accord with the current 
observational results, we have fl^F = 1 throughout the evolution, i.e. a scale-dependent 
Einstein-de Sitter model: 



3 a £ + i^>& = . («e)'_?^ = , (a . 6) 

a v 3a 6 V \a v J a v 



Suppose now that we 'fit' a standard Einstein-de Sitter model on some given scale D 
to observational data, we would be in the position to evaluate the physical 'parameters' 
on that scale to be = 1 — 4f2g, i.e. for ~ 1/3 today, we would conclude that 
there must be backreaction at work and it should be positive (negative kinematical 
backreaction mimicking a 'kinematical dark matter' source), fig ~ 1/6, and that the 
(physical) curvature parameter would be positive too (negative averaged curvature), 
fi^ = 3fig ~ 1/2. We emphasise that the Friedmannian curvature parameter is assumed 
to vanish, which demonstrates that it has nothing to do with the (evolving) averaged 
scalar curvature of the inhomogeneous model. 

The above procedure exemplifies the possibility of constructing a (non-naive) 
Friedmannian fitting model. It, however, assumes that, regionally, the model obeys an 
Einstein-de Sitter kinematics (unrelated to an underlying FRW metric), and r = 1/3 
is 'typical' for the regional Universe. Both are not in accord with what we expect. We 
would 'fit' a Friedmannian model with a cosmological constant, which is the currently 
held view of the 'concordance model', 

nZ F + + = i , (a.7) 

then we would have to superimpose a solution fig 2 with, e.g., r = — 1/3 to the above 
solution (this could still be interpreted as a deviation from a representative volume of a 
global model with r = 1/3). This implies, with fi^ F = Q% + 4fi|\ ~ 1/3 ; fi^ ~ 2/3, 
fi£ ~ 0, and 



that 



= -2VP Q2 ~ - . (A.9) 



2 
3 

We would find a positive backreaction with f2g 2 ~ — 1/3 (modeling now Dark Energy), 
and a negative averaged scalar curvature with ~ 1, indicating that the regional 
Universe should correspond to a 'void' within a global model with ~ 1/2. That our 
regional Universe could correspond to a regional 'void' has been discussed in a number 
of other papers, e.g. [73 EH1 EH , , EH, 0, EH- 

This example also demonstrates that we can construct cosmologies with different 
properties on global and regional scales by superimposing scaling solutions. 
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